UXarray for Advanced HEALPix Analysis & Visualization
In this section, you’ll learn:
Using the
uxarraypackage to perform advanced analysis operators over HEALPix data such as non-conservative zonal means, etc.
Related Documentation
Prerequisites
Concepts |
Importance |
Notes |
|---|---|---|
Necessary |
||
Necessary |
Time to learn: 30 minutes
import cartopy.crs as ccrs
import holoviews as hv
import intake
import uxarray as ux
Open data catalog
Tip
We assume, you have already gone over the previous section, UXarray for Basic HEALPix Statistics & Visualization. If not and if you need to learn about data catalogs in general and the data we will use throughout this notebook, we recommend to check that section first.:::
Let us open the online catalog from the WCRP’s Digital Earths Global Hackathon 2025 using intake and read the output of the ICON run ngc4008, which is stored in the HEALPix format:
# Final data catalog location (once hackathon website (https://digital-earths-global-hackathon.github.io/) updated)
# cat_url='https://digital-earths-global-hackathon.github.io/catalog/catalog.yaml'
# Interim data catalog location
cat_url = "https://raw.githubusercontent.com/digital-earths-global-hackathon/catalog/refs/heads/ncar/online/main.yaml"
cat = intake.open_catalog(cat_url)
model_run = cat.icon_ngc4008
We can look into a fine resolution dataset at zoome level = 10 in it as Xarray.Dataset:
ds = model_run(zoom=9, time="P1D").to_dask()
ds
/home/runner/miniconda3/envs/healpix-cookbook-dev/lib/python3.10/site-packages/intake_xarray/base.py:21: FutureWarning: The return type of `Dataset.dims` will be changed to return a set of dimension names in future, in order to be more consistent with `DataArray.dims`. To access a mapping from dimension names to lengths, please use `Dataset.sizes`.
'dims': dict(self._ds.dims),
<xarray.Dataset> Size: 232TB
Dimensions: (time: 10958, depth_half: 73,
cell: 3145728, level_full: 90,
crs: 1, depth_full: 72,
soil_depth_water_level: 5,
level_half: 91,
soil_depth_energy_level: 5)
Coordinates:
* crs (crs) float32 4B nan
* depth_full (depth_full) float32 288B 1.0 ... 5....
* depth_half (depth_half) float32 292B 0.0 ... 5....
* level_full (level_full) int32 360B 1 2 3 ... 89 90
* level_half (level_half) int32 364B 1 2 3 ... 90 91
* soil_depth_energy_level (soil_depth_energy_level) float32 20B ...
* soil_depth_water_level (soil_depth_water_level) float32 20B ...
* time (time) datetime64[ns] 88kB 2020-01-0...
Dimensions without coordinates: cell
Data variables: (12/103)
A_tracer_v_to (time, depth_half, cell) float32 10TB ...
FrshFlux_IceSalt (time, cell) float32 138GB ...
FrshFlux_TotalIce (time, cell) float32 138GB ...
Qbot (time, cell) float32 138GB ...
Qtop (time, cell) float32 138GB ...
Wind_Speed_10m (time, cell) float32 138GB ...
... ...
vas (time, cell) float32 138GB ...
w (time, depth_half, cell) float32 10TB ...
wa_phy (time, level_half, cell) float32 13TB ...
zg (level_full, cell) float32 1GB ...
zghalf (level_half, cell) float32 1GB ...
zos (time, cell) float32 138GB ...- time: 10958
- depth_half: 73
- cell: 3145728
- level_full: 90
- crs: 1
- depth_full: 72
- soil_depth_water_level: 5
- level_half: 91
- soil_depth_energy_level: 5
- crs(crs)float32nan
- grid_mapping_name :
- healpix
- healpix_nside :
- 512
- healpix_order :
- nest
array([nan], dtype=float32)
- depth_full(depth_full)float321.0 3.1 ... 5.546e+03 5.816e+03
- axis :
- Z
- long_name :
- depth_below_sea
- positive :
- down
- standard_name :
- depth
- units :
- m
array([1.00000e+00, 3.10000e+00, 5.45000e+00, 8.10000e+00, 1.10500e+01, 1.43500e+01, 1.80500e+01, 2.22000e+01, 2.68500e+01, 3.20000e+01, 3.76500e+01, 4.38000e+01, 5.05500e+01, 5.79500e+01, 6.60000e+01, 7.48000e+01, 8.44500e+01, 9.50000e+01, 1.06500e+02, 1.19100e+02, 1.32900e+02, 1.47950e+02, 1.64350e+02, 1.82250e+02, 2.01800e+02, 2.23150e+02, 2.46450e+02, 2.71850e+02, 2.99550e+02, 3.29750e+02, 3.62650e+02, 3.98450e+02, 4.37400e+02, 4.79250e+02, 5.23450e+02, 5.69500e+02, 6.17100e+02, 6.66300e+02, 7.17150e+02, 7.69700e+02, 8.24000e+02, 8.80100e+02, 9.38050e+02, 9.97900e+02, 1.05975e+03, 1.12450e+03, 1.19410e+03, 1.27070e+03, 1.35630e+03, 1.45170e+03, 1.55680e+03, 1.67195e+03, 1.79645e+03, 1.92950e+03, 2.07140e+03, 2.22245e+03, 2.38290e+03, 2.55290e+03, 2.73250e+03, 2.92175e+03, 3.12060e+03, 3.32885e+03, 3.54625e+03, 3.77245e+03, 4.00695e+03, 4.24915e+03, 4.49830e+03, 4.75355e+03, 5.01400e+03, 5.27860e+03, 5.54625e+03, 5.81575e+03], dtype=float32) - depth_half(depth_half)float320.0 2.0 4.2 ... 5.681e+03 5.951e+03
- axis :
- Z
- long_name :
- depth_below_sea
- positive :
- down
- standard_name :
- depth
- units :
- m
array([0.0000e+00, 2.0000e+00, 4.2000e+00, 6.7000e+00, 9.5000e+00, 1.2600e+01, 1.6100e+01, 2.0000e+01, 2.4400e+01, 2.9300e+01, 3.4700e+01, 4.0600e+01, 4.7000e+01, 5.4100e+01, 6.1800e+01, 7.0200e+01, 7.9400e+01, 8.9500e+01, 1.0050e+02, 1.1250e+02, 1.2570e+02, 1.4010e+02, 1.5580e+02, 1.7290e+02, 1.9160e+02, 2.1200e+02, 2.3430e+02, 2.5860e+02, 2.8510e+02, 3.1400e+02, 3.4550e+02, 3.7980e+02, 4.1710e+02, 4.5770e+02, 5.0080e+02, 5.4610e+02, 5.9290e+02, 6.4130e+02, 6.9130e+02, 7.4300e+02, 7.9640e+02, 8.5160e+02, 9.0860e+02, 9.6750e+02, 1.0283e+03, 1.0912e+03, 1.1578e+03, 1.2304e+03, 1.3110e+03, 1.4016e+03, 1.5018e+03, 1.6118e+03, 1.7321e+03, 1.8608e+03, 1.9982e+03, 2.1446e+03, 2.3003e+03, 2.4655e+03, 2.6403e+03, 2.8247e+03, 3.0188e+03, 3.2224e+03, 3.4353e+03, 3.6572e+03, 3.8877e+03, 4.1262e+03, 4.3721e+03, 4.6245e+03, 4.8826e+03, 5.1454e+03, 5.4118e+03, 5.6807e+03, 5.9508e+03], dtype=float32) - level_full(level_full)int321 2 3 4 5 6 7 ... 85 86 87 88 89 90
- axis :
- Z
- long_name :
- generalized_height
- positive :
- down
- standard_name :
- height
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90], dtype=int32) - level_half(level_half)int321 2 3 4 5 6 7 ... 86 87 88 89 90 91
- axis :
- Z
- long_name :
- generalized_height
- positive :
- down
- standard_name :
- height
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91], dtype=int32) - soil_depth_energy_level(soil_depth_energy_level)float320.0325 0.192 0.7755 2.683 6.984
- axis :
- Z
- long_name :
- depth_below_land
- positive :
- down
- units :
- m
array([0.0325, 0.192 , 0.7755, 2.683 , 6.984 ], dtype=float32)
- soil_depth_water_level(soil_depth_water_level)float320.0325 0.192 0.7755 2.683 6.984
- axis :
- Z
- long_name :
- depth_below_land
- positive :
- down
- units :
- m
array([0.0325, 0.192 , 0.7755, 2.683 , 6.984 ], dtype=float32)
- time(time)datetime64[ns]2020-01-02 ... 2050-01-01
- axis :
- T
array(['2020-01-02T00:00:00.000000000', '2020-01-03T00:00:00.000000000', '2020-01-04T00:00:00.000000000', ..., '2049-12-30T00:00:00.000000000', '2049-12-31T00:00:00.000000000', '2050-01-01T00:00:00.000000000'], shape=(10958,), dtype='datetime64[ns]')
- A_tracer_v_to(time, depth_half, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- sea water vertical diffusivity
- standard_name :
- ocean_vertical_diffusivity
- units :
- m2 s-1
- vgrid :
- depth_below_sea_half
[2516374781952 values with dtype=float32]
- FrshFlux_IceSalt(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Salt volume flux due to sea ice change
- units :
- psu m s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- FrshFlux_TotalIce(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Freshwater Flux due to Sea Ice Change
- units :
- m s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- Qbot(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Conductive heat flux at ice-ocean interface
- units :
- W/m^2
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- Qtop(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Energy flux available for surface melting
- units :
- W/m^2
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- Wind_Speed_10m(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Wind Speed at 10m height
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_FrshFlux_Evaporation(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_FrshFlux_Precipitation(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_FrshFlux_Runoff(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_FrshFlux_SnowFall(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_Latent(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_LongWave(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_Sensible(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_ShortWave(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_Total(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_stress_x(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_stress_xw(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_stress_y(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_stress_yw(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- cli(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific cloud ice content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- clivi(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated cloud ice
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- cllvi(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated cloud water
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- clw(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific cloud water content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- conc(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- ice concentration in each ice class
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- delhi(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Change in ice mean thickness due to thermodynamic effects
- units :
- m
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- dzghalf(level_full, cell)float32...
- long_name :
- vertical extent of model full layer
- standard_name :
- cell_thickness
- units :
- m
[283115520 values with dtype=float32]
- heatOceI(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Heat flux to ocean from the ice growth
- units :
- W/m^2
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- heatOceW(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Heat flux to ocean from the atmosphere
- units :
- W/m^2
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- heat_content_seaice(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- heat_content_snow(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- heat_content_total(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hfls(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- latent heat flux
- standard_name :
- surface_downward_latent_heat_flux
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hfss(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- sensible heat flux
- standard_name :
- surface_downward_sensible_heat_flux
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hi(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- ice thickness
- units :
- m
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- hs(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- snow thickness
- units :
- m
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- hus(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- Specific humidity
- standard_name :
- specific_humidity
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- hydro_canopy_cond_limited_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_discharge_ocean_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- units :
- m3 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_drainage_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- units :
- kg m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_ice_soil_sl_box(time, soil_depth_water_level, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Ice content in soil layers
- units :
- m
- vgrid :
- soil_depth_water
[172354437120 values with dtype=float32]
- hydro_runoff_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- surface runoff
- units :
- kg m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_snow_soil_dens_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Density of snow on soil
- units :
- kg m-3
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_transpiration_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Transpiration from surface
- units :
- kg m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_weq_snow_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Water content of snow reservoir on surface
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_wtr_soil_sl_box(time, soil_depth_water_level, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Water content in soil layers
- units :
- m
- vgrid :
- soil_depth_water
[172354437120 values with dtype=float32]
- ice_u(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- zonal velocity
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- ice_v(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- meridional velocity
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- mlotst(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- standard_name :
- ocean_mixed_layer_thickness_defined_by_sigma_t
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- newice(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- new ice growth in open water
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- ocean_fraction_depth_full(depth_full, cell)float32...
- long_name :
- fraction of cell covered by ocean
- standard_name :
- ocean_area_fraction
[226492416 values with dtype=float32]
- ocean_fraction_depth_half(depth_half, cell)float32...
- long_name :
- fraction of cell covered by ocean
- standard_name :
- ocean_area_fraction
[229638144 values with dtype=float32]
- ocean_fraction_surface(cell)float32...
- long_name :
- fraction of cell covered by ocean
- standard_name :
- ocean_area_fraction
[3145728 values with dtype=float32]
- pfull(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- standard_name :
- air_pressure
- units :
- Pa
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- phalf(time, level_half, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- standard_name :
- air_pressure
- units :
- Pa
- vgrid :
- reference_half
[3136850755584 values with dtype=float32]
- pr(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- precipitation flux
- units :
- kg m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- pres_msl(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- mean sea level pressure
- units :
- Pa
- vgrid :
- meansea
[34470887424 values with dtype=float32]
- pres_sfc(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface pressure
- standard_name :
- surface_air_pressure
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- prls(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- prw(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated water vapour
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- qg(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific graupel content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- qgvi(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated graupel
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- qr(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific rain content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- qrvi(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated rain
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- qs(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific snow content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- qsvi(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated snow
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- qv2m(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific humidity in 2m
- units :
- kg kg-1
- vgrid :
- height_2m
[34470887424 values with dtype=float32]
- rlds(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface downwelling longwave radiation
- standard_name :
- surface_downwelling_longwave_flux_in_air
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- rlus(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface upwelling longwave radiation
- standard_name :
- surface_upwelling_longwave_flux_in_air
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- rlut(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- toa outgoing longwave radiation
- standard_name :
- toa_outgoing_longwave_flux
- units :
- W m-2
- vgrid :
- toa
[34470887424 values with dtype=float32]
- rsds(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface downwelling shortwave radiation
- standard_name :
- surface_downwelling_shortwave_flux_in_air
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- rsdt(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- toa incident shortwave radiation
- standard_name :
- toa_incoming_shortwave_flux
- units :
- W m-2
- vgrid :
- toa
[34470887424 values with dtype=float32]
- rsus(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface upwelling shortwave radiation
- standard_name :
- surface_upwelling_shortwave_flux_in_air
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- rsut(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- toa outgoing shortwave radiation
- standard_name :
- toa_outgoing_shortwave_flux
- units :
- W m-2
- vgrid :
- toa
[34470887424 values with dtype=float32]
- seb_forc_hflx_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- additional heat flux at surface
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- seb_heat_cap_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- surface layer heat capacity
- units :
- J m-2 K-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- sfcwind(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- 10m windspeed
- units :
- m s-1
- vgrid :
- height_10m
[34470887424 values with dtype=float32]
- sic(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- fraction of ocean covered by sea ice
- vgrid :
- surface
[34470887424 values with dtype=float32]
- sit(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- sea ice thickness
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- snow_to_ice(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- amount of snow that is transformed to ice
- units :
- m
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- so(time, depth_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- sea water salinity
- standard_name :
- sea_water_salinity
- units :
- psu
- vgrid :
- depth_below_sea
[2481903894528 values with dtype=float32]
- sse_grnd_hflx_old_box(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Ground heat flux (old)
- units :
- J m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- sse_t_soil_sl_box(time, soil_depth_energy_level, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- standard_name :
- soil_temperature
- units :
- K
- vgrid :
- soil_depth_energy
[172354437120 values with dtype=float32]
- stretch_c(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- zstar surface stretch at cell center
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- ta(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- Temperature
- standard_name :
- air_temperature
- units :
- K
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- tas(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- temperature in 2m
- standard_name :
- air_temperature
- units :
- K
- vgrid :
- height_2m
[34470887424 values with dtype=float32]
- tauu(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- u-momentum flux at the surface
- units :
- N m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- tauv(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- v-momentum flux at the surface
- units :
- N m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- tke(time, depth_half, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- turbulent kinetic energy
- units :
- m2 s-2
- vgrid :
- depth_below_sea_half
[2516374781952 values with dtype=float32]
- to(time, depth_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- sea water potential temperature
- standard_name :
- sea_water_potential_temperature
- units :
- degC
- vgrid :
- depth_below_sea
[2481903894528 values with dtype=float32]
- ts(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface temperature
- standard_name :
- surface_temperature
- units :
- K
- vgrid :
- surface
[34470887424 values with dtype=float32]
- u(time, depth_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- u zonal velocity component
- units :
- m/s
- vgrid :
- depth_below_sea
[2481903894528 values with dtype=float32]
- u_vint(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- vertically integrated barotropic zonal velocity
- units :
- m2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- ua(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- Zonal wind
- standard_name :
- eastward_wind
- units :
- m s-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- uas(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- zonal wind in 10m
- units :
- m s-1
- vgrid :
- height_10m
[34470887424 values with dtype=float32]
- v(time, depth_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- v meridional velocity component
- units :
- m/s
- vgrid :
- depth_below_sea
[2481903894528 values with dtype=float32]
- v_vint(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- vertically integrated meridional zonal velocity
- units :
- m2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- va(time, level_full, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- Meridional wind
- standard_name :
- northward_wind
- units :
- m s-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- vas(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- meridional wind in 10m
- units :
- m s-1
- vgrid :
- height_10m
[34470887424 values with dtype=float32]
- w(time, depth_half, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- vertical velocity at cells
- units :
- m/s
- vgrid :
- depth_below_sea_half
[2516374781952 values with dtype=float32]
- wa_phy(time, level_half, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertical velocity
- units :
- m s-1
- vgrid :
- reference_half
[3136850755584 values with dtype=float32]
- zg(level_full, cell)float32...
- long_name :
- geometric height at full level center
- standard_name :
- height
- units :
- m
[283115520 values with dtype=float32]
- zghalf(level_half, cell)float32...
- long_name :
- geometric height at half level center
- standard_name :
- height
- units :
- m
[286261248 values with dtype=float32]
- zos(time, cell)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- zstar sfc elevation at cell center
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- crsPandasIndex
PandasIndex(Index([nan], dtype='float32', name='crs'))
- depth_fullPandasIndex
PandasIndex(Index([ 1.0, 3.0999999046325684, 5.449999809265137, 8.100000381469727, 11.050000190734863, 14.350000381469727, 18.049999237060547, 22.200000762939453, 26.850000381469727, 32.0, 37.650001525878906, 43.79999923706055, 50.54999923706055, 57.95000076293945, 66.0, 74.80000305175781, 84.44999694824219, 95.0, 106.5, 119.0999984741211, 132.89999389648438, 147.9499969482422, 164.35000610351562, 182.25, 201.8000030517578, 223.14999389648438, 246.4499969482422, 271.8500061035156, 299.54998779296875, 329.75, 362.6499938964844, 398.45001220703125, 437.3999938964844, 479.25, 523.4500122070312, 569.5, 617.0999755859375, 666.2999877929688, 717.1500244140625, 769.7000122070312, 824.0, 880.0999755859375, 938.0499877929688, 997.9000244140625, 1059.75, 1124.5, 1194.0999755859375, 1270.699951171875, 1356.300048828125, 1451.699951171875, 1556.800048828125, 1671.949951171875, 1796.449951171875, 1929.5, 2071.39990234375, 2222.449951171875, 2382.89990234375, 2552.89990234375, 2732.5, 2921.75, 3120.60009765625, 3328.85009765625, 3546.25, 3772.449951171875, 4006.949951171875, 4249.14990234375, 4498.2998046875, 4753.5498046875, 5014.0, 5278.60009765625, 5546.25, 5815.75], dtype='float32', name='depth_full')) - depth_halfPandasIndex
PandasIndex(Index([ 0.0, 2.0, 4.199999809265137, 6.699999809265137, 9.5, 12.600000381469727, 16.100000381469727, 20.0, 24.399999618530273, 29.299999237060547, 34.70000076293945, 40.599998474121094, 47.0, 54.099998474121094, 61.79999923706055, 70.19999694824219, 79.4000015258789, 89.5, 100.5, 112.5, 125.69999694824219, 140.10000610351562, 155.8000030517578, 172.89999389648438, 191.60000610351562, 212.0, 234.3000030517578, 258.6000061035156, 285.1000061035156, 314.0, 345.5, 379.79998779296875, 417.1000061035156, 457.70001220703125, 500.79998779296875, 546.0999755859375, 592.9000244140625, 641.2999877929688, 691.2999877929688, 743.0, 796.4000244140625, 851.5999755859375, 908.5999755859375, 967.5, 1028.300048828125, 1091.199951171875, 1157.800048828125, 1230.4000244140625, 1311.0, 1401.5999755859375, 1501.800048828125, 1611.800048828125, 1732.0999755859375, 1860.800048828125, 1998.199951171875, 2144.60009765625, 2300.300048828125, 2465.5, 2640.300048828125, 2824.699951171875, 3018.800048828125, 3222.39990234375, 3435.300048828125, 3657.199951171875, 3887.699951171875, 4126.2001953125, 4372.10009765625, 4624.5, 4882.60009765625, 5145.39990234375, 5411.7998046875, 5680.7001953125, 5950.7998046875], dtype='float32', name='depth_half')) - level_fullPandasIndex
PandasIndex(Index([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90], dtype='int32', name='level_full')) - level_halfPandasIndex
PandasIndex(Index([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91], dtype='int32', name='level_half')) - soil_depth_energy_levelPandasIndex
PandasIndex(Index([0.032499998807907104, 0.19200000166893005, 0.7754999995231628, 2.683000087738037, 6.984000205993652], dtype='float32', name='soil_depth_energy_level')) - soil_depth_water_levelPandasIndex
PandasIndex(Index([0.032499998807907104, 0.19200000166893005, 0.7754999995231628, 2.683000087738037, 6.984000205993652], dtype='float32', name='soil_depth_water_level')) - timePandasIndex
PandasIndex(DatetimeIndex(['2020-01-02', '2020-01-03', '2020-01-04', '2020-01-05', '2020-01-06', '2020-01-07', '2020-01-08', '2020-01-09', '2020-01-10', '2020-01-11', ... '2049-12-23', '2049-12-24', '2049-12-25', '2049-12-26', '2049-12-27', '2049-12-28', '2049-12-29', '2049-12-30', '2049-12-31', '2050-01-01'], dtype='datetime64[ns]', name='time', length=10958, freq=None))
Create UXarray Datasets from HEALPix
We can use from_healpix as follows to open a HEALPix grid from xarray.Dataset:
uxds = ux.UxDataset.from_healpix(ds)
uxds
<xarray.UxDataset> Size: 232TB
Dimensions: (time: 10958, depth_half: 73,
n_face: 3145728, level_full: 90,
crs: 1, depth_full: 72,
soil_depth_water_level: 5,
level_half: 91,
soil_depth_energy_level: 5)
Coordinates:
* crs (crs) float32 4B nan
* depth_full (depth_full) float32 288B 1.0 ... 5....
* depth_half (depth_half) float32 292B 0.0 ... 5....
* level_full (level_full) int32 360B 1 2 3 ... 89 90
* level_half (level_half) int32 364B 1 2 3 ... 90 91
* soil_depth_energy_level (soil_depth_energy_level) float32 20B ...
* soil_depth_water_level (soil_depth_water_level) float32 20B ...
* time (time) datetime64[ns] 88kB 2020-01-0...
Dimensions without coordinates: n_face
Data variables: (12/103)
A_tracer_v_to (time, depth_half, n_face) float32 10TB ...
FrshFlux_IceSalt (time, n_face) float32 138GB ...
FrshFlux_TotalIce (time, n_face) float32 138GB ...
Qbot (time, n_face) float32 138GB ...
Qtop (time, n_face) float32 138GB ...
Wind_Speed_10m (time, n_face) float32 138GB ...
... ...
vas (time, n_face) float32 138GB ...
w (time, depth_half, n_face) float32 10TB ...
wa_phy (time, level_half, n_face) float32 13TB ...
zg (level_full, n_face) float32 1GB ...
zghalf (level_half, n_face) float32 1GB ...
zos (time, n_face) float32 138GB ...<xarray.UxDataset> Size: 232TB
Dimensions: (time: 10958, depth_half: 73,
n_face: 3145728, level_full: 90,
crs: 1, depth_full: 72,
soil_depth_water_level: 5,
level_half: 91,
soil_depth_energy_level: 5)
Coordinates:
* crs (crs) float32 4B nan
* depth_full (depth_full) float32 288B 1.0 ... 5....
* depth_half (depth_half) float32 292B 0.0 ... 5....
* level_full (level_full) int32 360B 1 2 3 ... 89 90
* level_half (level_half) int32 364B 1 2 3 ... 90 91
* soil_depth_energy_level (soil_depth_energy_level) float32 20B ...
* soil_depth_water_level (soil_depth_water_level) float32 20B ...
* time (time) datetime64[ns] 88kB 2020-01-0...
Dimensions without coordinates: n_face
Data variables: (12/103)
A_tracer_v_to (time, depth_half, n_face) float32 10TB ...
FrshFlux_IceSalt (time, n_face) float32 138GB ...
FrshFlux_TotalIce (time, n_face) float32 138GB ...
Qbot (time, n_face) float32 138GB ...
Qtop (time, n_face) float32 138GB ...
Wind_Speed_10m (time, n_face) float32 138GB ...
... ...
vas (time, n_face) float32 138GB ...
w (time, depth_half, n_face) float32 10TB ...
wa_phy (time, level_half, n_face) float32 13TB ...
zg (level_full, n_face) float32 1GB ...
zghalf (level_half, n_face) float32 1GB ...
zos (time, n_face) float32 138GB ...- time: 10958
- depth_half: 73
- n_face: 3145728
- level_full: 90
- crs: 1
- depth_full: 72
- soil_depth_water_level: 5
- level_half: 91
- soil_depth_energy_level: 5
- crs(crs)float32nan
- grid_mapping_name :
- healpix
- healpix_nside :
- 512
- healpix_order :
- nest
array([nan], dtype=float32)
- depth_full(depth_full)float321.0 3.1 ... 5.546e+03 5.816e+03
- axis :
- Z
- long_name :
- depth_below_sea
- positive :
- down
- standard_name :
- depth
- units :
- m
array([1.00000e+00, 3.10000e+00, 5.45000e+00, 8.10000e+00, 1.10500e+01, 1.43500e+01, 1.80500e+01, 2.22000e+01, 2.68500e+01, 3.20000e+01, 3.76500e+01, 4.38000e+01, 5.05500e+01, 5.79500e+01, 6.60000e+01, 7.48000e+01, 8.44500e+01, 9.50000e+01, 1.06500e+02, 1.19100e+02, 1.32900e+02, 1.47950e+02, 1.64350e+02, 1.82250e+02, 2.01800e+02, 2.23150e+02, 2.46450e+02, 2.71850e+02, 2.99550e+02, 3.29750e+02, 3.62650e+02, 3.98450e+02, 4.37400e+02, 4.79250e+02, 5.23450e+02, 5.69500e+02, 6.17100e+02, 6.66300e+02, 7.17150e+02, 7.69700e+02, 8.24000e+02, 8.80100e+02, 9.38050e+02, 9.97900e+02, 1.05975e+03, 1.12450e+03, 1.19410e+03, 1.27070e+03, 1.35630e+03, 1.45170e+03, 1.55680e+03, 1.67195e+03, 1.79645e+03, 1.92950e+03, 2.07140e+03, 2.22245e+03, 2.38290e+03, 2.55290e+03, 2.73250e+03, 2.92175e+03, 3.12060e+03, 3.32885e+03, 3.54625e+03, 3.77245e+03, 4.00695e+03, 4.24915e+03, 4.49830e+03, 4.75355e+03, 5.01400e+03, 5.27860e+03, 5.54625e+03, 5.81575e+03], dtype=float32) - depth_half(depth_half)float320.0 2.0 4.2 ... 5.681e+03 5.951e+03
- axis :
- Z
- long_name :
- depth_below_sea
- positive :
- down
- standard_name :
- depth
- units :
- m
array([0.0000e+00, 2.0000e+00, 4.2000e+00, 6.7000e+00, 9.5000e+00, 1.2600e+01, 1.6100e+01, 2.0000e+01, 2.4400e+01, 2.9300e+01, 3.4700e+01, 4.0600e+01, 4.7000e+01, 5.4100e+01, 6.1800e+01, 7.0200e+01, 7.9400e+01, 8.9500e+01, 1.0050e+02, 1.1250e+02, 1.2570e+02, 1.4010e+02, 1.5580e+02, 1.7290e+02, 1.9160e+02, 2.1200e+02, 2.3430e+02, 2.5860e+02, 2.8510e+02, 3.1400e+02, 3.4550e+02, 3.7980e+02, 4.1710e+02, 4.5770e+02, 5.0080e+02, 5.4610e+02, 5.9290e+02, 6.4130e+02, 6.9130e+02, 7.4300e+02, 7.9640e+02, 8.5160e+02, 9.0860e+02, 9.6750e+02, 1.0283e+03, 1.0912e+03, 1.1578e+03, 1.2304e+03, 1.3110e+03, 1.4016e+03, 1.5018e+03, 1.6118e+03, 1.7321e+03, 1.8608e+03, 1.9982e+03, 2.1446e+03, 2.3003e+03, 2.4655e+03, 2.6403e+03, 2.8247e+03, 3.0188e+03, 3.2224e+03, 3.4353e+03, 3.6572e+03, 3.8877e+03, 4.1262e+03, 4.3721e+03, 4.6245e+03, 4.8826e+03, 5.1454e+03, 5.4118e+03, 5.6807e+03, 5.9508e+03], dtype=float32) - level_full(level_full)int321 2 3 4 5 6 7 ... 85 86 87 88 89 90
- axis :
- Z
- long_name :
- generalized_height
- positive :
- down
- standard_name :
- height
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90], dtype=int32) - level_half(level_half)int321 2 3 4 5 6 7 ... 86 87 88 89 90 91
- axis :
- Z
- long_name :
- generalized_height
- positive :
- down
- standard_name :
- height
array([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91], dtype=int32) - soil_depth_energy_level(soil_depth_energy_level)float320.0325 0.192 0.7755 2.683 6.984
- axis :
- Z
- long_name :
- depth_below_land
- positive :
- down
- units :
- m
array([0.0325, 0.192 , 0.7755, 2.683 , 6.984 ], dtype=float32)
- soil_depth_water_level(soil_depth_water_level)float320.0325 0.192 0.7755 2.683 6.984
- axis :
- Z
- long_name :
- depth_below_land
- positive :
- down
- units :
- m
array([0.0325, 0.192 , 0.7755, 2.683 , 6.984 ], dtype=float32)
- time(time)datetime64[ns]2020-01-02 ... 2050-01-01
- axis :
- T
array(['2020-01-02T00:00:00.000000000', '2020-01-03T00:00:00.000000000', '2020-01-04T00:00:00.000000000', ..., '2049-12-30T00:00:00.000000000', '2049-12-31T00:00:00.000000000', '2050-01-01T00:00:00.000000000'], shape=(10958,), dtype='datetime64[ns]')
- A_tracer_v_to(time, depth_half, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- sea water vertical diffusivity
- standard_name :
- ocean_vertical_diffusivity
- units :
- m2 s-1
- vgrid :
- depth_below_sea_half
[2516374781952 values with dtype=float32]
- FrshFlux_IceSalt(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Salt volume flux due to sea ice change
- units :
- psu m s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- FrshFlux_TotalIce(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Freshwater Flux due to Sea Ice Change
- units :
- m s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- Qbot(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Conductive heat flux at ice-ocean interface
- units :
- W/m^2
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- Qtop(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Energy flux available for surface melting
- units :
- W/m^2
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- Wind_Speed_10m(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Wind Speed at 10m height
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_FrshFlux_Evaporation(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_FrshFlux_Precipitation(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_FrshFlux_Runoff(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_FrshFlux_SnowFall(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_Latent(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_LongWave(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_Sensible(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_ShortWave(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_HeatFlux_Total(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- W/m2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_stress_x(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_stress_xw(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_stress_y(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- atmos_fluxes_stress_yw(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- cli(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific cloud ice content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- clivi(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated cloud ice
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- cllvi(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated cloud water
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- clw(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific cloud water content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- conc(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- ice concentration in each ice class
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- delhi(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Change in ice mean thickness due to thermodynamic effects
- units :
- m
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- dzghalf(level_full, n_face)float32...
- long_name :
- vertical extent of model full layer
- standard_name :
- cell_thickness
- units :
- m
[283115520 values with dtype=float32]
- heatOceI(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Heat flux to ocean from the ice growth
- units :
- W/m^2
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- heatOceW(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- Heat flux to ocean from the atmosphere
- units :
- W/m^2
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- heat_content_seaice(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- heat_content_snow(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- heat_content_total(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hfls(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- latent heat flux
- standard_name :
- surface_downward_latent_heat_flux
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hfss(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- sensible heat flux
- standard_name :
- surface_downward_sensible_heat_flux
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hi(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- ice thickness
- units :
- m
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- hs(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- snow thickness
- units :
- m
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- hus(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- Specific humidity
- standard_name :
- specific_humidity
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- hydro_canopy_cond_limited_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_discharge_ocean_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- units :
- m3 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_drainage_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- units :
- kg m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_ice_soil_sl_box(time, soil_depth_water_level, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Ice content in soil layers
- units :
- m
- vgrid :
- soil_depth_water
[172354437120 values with dtype=float32]
- hydro_runoff_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- surface runoff
- units :
- kg m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_snow_soil_dens_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Density of snow on soil
- units :
- kg m-3
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_transpiration_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Transpiration from surface
- units :
- kg m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_weq_snow_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Water content of snow reservoir on surface
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- hydro_wtr_soil_sl_box(time, soil_depth_water_level, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Water content in soil layers
- units :
- m
- vgrid :
- soil_depth_water
[172354437120 values with dtype=float32]
- ice_u(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- zonal velocity
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- ice_v(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- meridional velocity
- units :
- m/s
- vgrid :
- surface
[34470887424 values with dtype=float32]
- mlotst(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- standard_name :
- ocean_mixed_layer_thickness_defined_by_sigma_t
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- newice(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- new ice growth in open water
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- ocean_fraction_depth_full(depth_full, n_face)float32...
- long_name :
- fraction of cell covered by ocean
- standard_name :
- ocean_area_fraction
[226492416 values with dtype=float32]
- ocean_fraction_depth_half(depth_half, n_face)float32...
- long_name :
- fraction of cell covered by ocean
- standard_name :
- ocean_area_fraction
[229638144 values with dtype=float32]
- ocean_fraction_surface(n_face)float32...
- long_name :
- fraction of cell covered by ocean
- standard_name :
- ocean_area_fraction
[3145728 values with dtype=float32]
- pfull(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- standard_name :
- air_pressure
- units :
- Pa
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- phalf(time, level_half, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- standard_name :
- air_pressure
- units :
- Pa
- vgrid :
- reference_half
[3136850755584 values with dtype=float32]
- pr(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- precipitation flux
- units :
- kg m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- pres_msl(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- mean sea level pressure
- units :
- Pa
- vgrid :
- meansea
[34470887424 values with dtype=float32]
- pres_sfc(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface pressure
- standard_name :
- surface_air_pressure
- units :
- Pa
- vgrid :
- surface
[34470887424 values with dtype=float32]
- prls(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- vgrid :
- surface
[34470887424 values with dtype=float32]
- prw(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated water vapour
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- qg(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific graupel content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- qgvi(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated graupel
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- qr(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific rain content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- qrvi(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated rain
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- qs(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific snow content
- units :
- kg kg-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- qsvi(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertically integrated snow
- units :
- kg m-2
- vgrid :
- atmosphere
[34470887424 values with dtype=float32]
- qv2m(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- specific humidity in 2m
- units :
- kg kg-1
- vgrid :
- height_2m
[34470887424 values with dtype=float32]
- rlds(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface downwelling longwave radiation
- standard_name :
- surface_downwelling_longwave_flux_in_air
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- rlus(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface upwelling longwave radiation
- standard_name :
- surface_upwelling_longwave_flux_in_air
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- rlut(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- toa outgoing longwave radiation
- standard_name :
- toa_outgoing_longwave_flux
- units :
- W m-2
- vgrid :
- toa
[34470887424 values with dtype=float32]
- rsds(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface downwelling shortwave radiation
- standard_name :
- surface_downwelling_shortwave_flux_in_air
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- rsdt(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- toa incident shortwave radiation
- standard_name :
- toa_incoming_shortwave_flux
- units :
- W m-2
- vgrid :
- toa
[34470887424 values with dtype=float32]
- rsus(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface upwelling shortwave radiation
- standard_name :
- surface_upwelling_shortwave_flux_in_air
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- rsut(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- toa outgoing shortwave radiation
- standard_name :
- toa_outgoing_shortwave_flux
- units :
- W m-2
- vgrid :
- toa
[34470887424 values with dtype=float32]
- seb_forc_hflx_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- additional heat flux at surface
- units :
- W m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- seb_heat_cap_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- surface layer heat capacity
- units :
- J m-2 K-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- sfcwind(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- 10m windspeed
- units :
- m s-1
- vgrid :
- height_10m
[34470887424 values with dtype=float32]
- sic(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- fraction of ocean covered by sea ice
- vgrid :
- surface
[34470887424 values with dtype=float32]
- sit(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- sea ice thickness
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- snow_to_ice(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- amount of snow that is transformed to ice
- units :
- m
- vgrid :
- generic_ice
[34470887424 values with dtype=float32]
- so(time, depth_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- sea water salinity
- standard_name :
- sea_water_salinity
- units :
- psu
- vgrid :
- depth_below_sea
[2481903894528 values with dtype=float32]
- sse_grnd_hflx_old_box(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- long_name :
- Ground heat flux (old)
- units :
- J m-2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- sse_t_soil_sl_box(time, soil_depth_energy_level, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- jsbach
- grid_mapping :
- crs
- standard_name :
- soil_temperature
- units :
- K
- vgrid :
- soil_depth_energy
[172354437120 values with dtype=float32]
- stretch_c(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- zstar surface stretch at cell center
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- ta(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- Temperature
- standard_name :
- air_temperature
- units :
- K
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- tas(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- temperature in 2m
- standard_name :
- air_temperature
- units :
- K
- vgrid :
- height_2m
[34470887424 values with dtype=float32]
- tauu(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- u-momentum flux at the surface
- units :
- N m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- tauv(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- v-momentum flux at the surface
- units :
- N m-2
- vgrid :
- surface
[34470887424 values with dtype=float32]
- tke(time, depth_half, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- turbulent kinetic energy
- units :
- m2 s-2
- vgrid :
- depth_below_sea_half
[2516374781952 values with dtype=float32]
- to(time, depth_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- sea water potential temperature
- standard_name :
- sea_water_potential_temperature
- units :
- degC
- vgrid :
- depth_below_sea
[2481903894528 values with dtype=float32]
- ts(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- surface temperature
- standard_name :
- surface_temperature
- units :
- K
- vgrid :
- surface
[34470887424 values with dtype=float32]
- u(time, depth_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- u zonal velocity component
- units :
- m/s
- vgrid :
- depth_below_sea
[2481903894528 values with dtype=float32]
- u_vint(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- vertically integrated barotropic zonal velocity
- units :
- m2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- ua(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- Zonal wind
- standard_name :
- eastward_wind
- units :
- m s-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- uas(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- zonal wind in 10m
- units :
- m s-1
- vgrid :
- height_10m
[34470887424 values with dtype=float32]
- v(time, depth_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- v meridional velocity component
- units :
- m/s
- vgrid :
- depth_below_sea
[2481903894528 values with dtype=float32]
- v_vint(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- vertically integrated meridional zonal velocity
- units :
- m2 s-1
- vgrid :
- surface
[34470887424 values with dtype=float32]
- va(time, level_full, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- Meridional wind
- standard_name :
- northward_wind
- units :
- m s-1
- vgrid :
- reference
[3102379868160 values with dtype=float32]
- vas(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- meridional wind in 10m
- units :
- m s-1
- vgrid :
- height_10m
[34470887424 values with dtype=float32]
- w(time, depth_half, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- vertical velocity at cells
- units :
- m/s
- vgrid :
- depth_below_sea_half
[2516374781952 values with dtype=float32]
- wa_phy(time, level_half, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- vertical velocity
- units :
- m s-1
- vgrid :
- reference_half
[3136850755584 values with dtype=float32]
- zg(level_full, n_face)float32...
- long_name :
- geometric height at full level center
- standard_name :
- height
- units :
- m
[283115520 values with dtype=float32]
- zghalf(level_half, n_face)float32...
- long_name :
- geometric height at half level center
- standard_name :
- height
- units :
- m
[286261248 values with dtype=float32]
- zos(time, n_face)float32...
- cell_methods :
- time: mean cell: mean
- component :
- ocean
- grid_mapping :
- crs
- long_name :
- zstar sfc elevation at cell center
- units :
- m
- vgrid :
- surface
[34470887424 values with dtype=float32]
- crsPandasIndex
PandasIndex(Index([nan], dtype='float32', name='crs'))
- depth_fullPandasIndex
PandasIndex(Index([ 1.0, 3.0999999046325684, 5.449999809265137, 8.100000381469727, 11.050000190734863, 14.350000381469727, 18.049999237060547, 22.200000762939453, 26.850000381469727, 32.0, 37.650001525878906, 43.79999923706055, 50.54999923706055, 57.95000076293945, 66.0, 74.80000305175781, 84.44999694824219, 95.0, 106.5, 119.0999984741211, 132.89999389648438, 147.9499969482422, 164.35000610351562, 182.25, 201.8000030517578, 223.14999389648438, 246.4499969482422, 271.8500061035156, 299.54998779296875, 329.75, 362.6499938964844, 398.45001220703125, 437.3999938964844, 479.25, 523.4500122070312, 569.5, 617.0999755859375, 666.2999877929688, 717.1500244140625, 769.7000122070312, 824.0, 880.0999755859375, 938.0499877929688, 997.9000244140625, 1059.75, 1124.5, 1194.0999755859375, 1270.699951171875, 1356.300048828125, 1451.699951171875, 1556.800048828125, 1671.949951171875, 1796.449951171875, 1929.5, 2071.39990234375, 2222.449951171875, 2382.89990234375, 2552.89990234375, 2732.5, 2921.75, 3120.60009765625, 3328.85009765625, 3546.25, 3772.449951171875, 4006.949951171875, 4249.14990234375, 4498.2998046875, 4753.5498046875, 5014.0, 5278.60009765625, 5546.25, 5815.75], dtype='float32', name='depth_full')) - depth_halfPandasIndex
PandasIndex(Index([ 0.0, 2.0, 4.199999809265137, 6.699999809265137, 9.5, 12.600000381469727, 16.100000381469727, 20.0, 24.399999618530273, 29.299999237060547, 34.70000076293945, 40.599998474121094, 47.0, 54.099998474121094, 61.79999923706055, 70.19999694824219, 79.4000015258789, 89.5, 100.5, 112.5, 125.69999694824219, 140.10000610351562, 155.8000030517578, 172.89999389648438, 191.60000610351562, 212.0, 234.3000030517578, 258.6000061035156, 285.1000061035156, 314.0, 345.5, 379.79998779296875, 417.1000061035156, 457.70001220703125, 500.79998779296875, 546.0999755859375, 592.9000244140625, 641.2999877929688, 691.2999877929688, 743.0, 796.4000244140625, 851.5999755859375, 908.5999755859375, 967.5, 1028.300048828125, 1091.199951171875, 1157.800048828125, 1230.4000244140625, 1311.0, 1401.5999755859375, 1501.800048828125, 1611.800048828125, 1732.0999755859375, 1860.800048828125, 1998.199951171875, 2144.60009765625, 2300.300048828125, 2465.5, 2640.300048828125, 2824.699951171875, 3018.800048828125, 3222.39990234375, 3435.300048828125, 3657.199951171875, 3887.699951171875, 4126.2001953125, 4372.10009765625, 4624.5, 4882.60009765625, 5145.39990234375, 5411.7998046875, 5680.7001953125, 5950.7998046875], dtype='float32', name='depth_half')) - level_fullPandasIndex
PandasIndex(Index([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90], dtype='int32', name='level_full')) - level_halfPandasIndex
PandasIndex(Index([ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84, 85, 86, 87, 88, 89, 90, 91], dtype='int32', name='level_half')) - soil_depth_energy_levelPandasIndex
PandasIndex(Index([0.032499998807907104, 0.19200000166893005, 0.7754999995231628, 2.683000087738037, 6.984000205993652], dtype='float32', name='soil_depth_energy_level')) - soil_depth_water_levelPandasIndex
PandasIndex(Index([0.032499998807907104, 0.19200000166893005, 0.7754999995231628, 2.683000087738037, 6.984000205993652], dtype='float32', name='soil_depth_water_level')) - timePandasIndex
PandasIndex(DatetimeIndex(['2020-01-02', '2020-01-03', '2020-01-04', '2020-01-05', '2020-01-06', '2020-01-07', '2020-01-08', '2020-01-09', '2020-01-10', '2020-01-11', ... '2049-12-23', '2049-12-24', '2049-12-25', '2049-12-26', '2049-12-27', '2049-12-28', '2049-12-29', '2049-12-30', '2049-12-31', '2050-01-01'], dtype='datetime64[ns]', name='time', length=10958, freq=None))
Show Grid Information
<uxarray.Grid> Original Grid Type: HEALPix Grid Dimensions: * n_face: 3145728 Grid Coordinates (Spherical): * face_lon: (3145728,) * face_lat: (3145728,) Grid Coordinates (Cartesian): Grid Connectivity Variables: Grid Descriptor Variables:
- n_face: 3145728
- face_lon(n_face)float6445.0 45.09 44.91 ... -45.09 -45.0
- standard_name :
- longitude
- long name :
- Longitude of the center of each face
- units :
- degrees_east
array([ 45. , 45.08789062, 44.91210938, ..., -44.91210938, -45.08789062, -45. ], shape=(3145728,)) - face_lat(n_face)float640.0746 0.1492 ... -0.1492 -0.0746
- standard_name :
- latitude
- long name :
- Latitude of the center of each face
- units :
- degrees_north
array([ 0.0746039 , 0.14920793, 0.14920793, ..., -0.14920793, -0.14920793, -0.0746039 ], shape=(3145728,))
- zoom :
- 9
- n_side :
- 512
- n_pix :
- 3145728
- nest :
- True
- source_grid_spec :
- HEALPix
Data variable of interest
Then let us pick a variable from the dataset, which will give us an uxarray.UxDataArray:
uxda = uxds["tas"]
uxda
<xarray.UxDataArray 'tas' (time: 10958, n_face: 3145728)> Size: 138GB
[34470887424 values with dtype=float32]
Coordinates:
* time (time) datetime64[ns] 88kB 2020-01-02 2020-01-03 ... 2050-01-01
Dimensions without coordinates: n_face
Attributes:
cell_methods: time: mean cell: mean
component: atmo
grid_mapping: crs
long_name: temperature in 2m
standard_name: air_temperature
units: K
vgrid: height_2m<xarray.UxDataArray 'tas' (time: 10958, n_face: 3145728)> Size: 138GB
[34470887424 values with dtype=float32]
Coordinates:
* time (time) datetime64[ns] 88kB 2020-01-02 2020-01-03 ... 2050-01-01
Dimensions without coordinates: n_face
Attributes:
cell_methods: time: mean cell: mean
component: atmo
grid_mapping: crs
long_name: temperature in 2m
standard_name: air_temperature
units: K
vgrid: height_2m- time: 10958
- n_face: 3145728
- ...
[34470887424 values with dtype=float32]
- time(time)datetime64[ns]2020-01-02 ... 2050-01-01
- axis :
- T
array(['2020-01-02T00:00:00.000000000', '2020-01-03T00:00:00.000000000', '2020-01-04T00:00:00.000000000', ..., '2049-12-30T00:00:00.000000000', '2049-12-31T00:00:00.000000000', '2050-01-01T00:00:00.000000000'], shape=(10958,), dtype='datetime64[ns]')
- timePandasIndex
PandasIndex(DatetimeIndex(['2020-01-02', '2020-01-03', '2020-01-04', '2020-01-05', '2020-01-06', '2020-01-07', '2020-01-08', '2020-01-09', '2020-01-10', '2020-01-11', ... '2049-12-23', '2049-12-24', '2049-12-25', '2049-12-26', '2049-12-27', '2049-12-28', '2049-12-29', '2049-12-30', '2049-12-31', '2050-01-01'], dtype='datetime64[ns]', name='time', length=10958, freq=None))
- cell_methods :
- time: mean cell: mean
- component :
- atmo
- grid_mapping :
- crs
- long_name :
- temperature in 2m
- standard_name :
- air_temperature
- units :
- K
- vgrid :
- height_2m
Show Grid Information
<uxarray.Grid> Original Grid Type: HEALPix Grid Dimensions: * n_face: 3145728 Grid Coordinates (Spherical): * face_lon: (3145728,) * face_lat: (3145728,) Grid Coordinates (Cartesian): Grid Connectivity Variables: Grid Descriptor Variables:
- n_face: 3145728
- face_lon(n_face)float6445.0 45.09 44.91 ... -45.09 -45.0
- standard_name :
- longitude
- long name :
- Longitude of the center of each face
- units :
- degrees_east
array([ 45. , 45.08789062, 44.91210938, ..., -44.91210938, -45.08789062, -45. ], shape=(3145728,)) - face_lat(n_face)float640.0746 0.1492 ... -0.1492 -0.0746
- standard_name :
- latitude
- long name :
- Latitude of the center of each face
- units :
- degrees_north
array([ 0.0746039 , 0.14920793, 0.14920793, ..., -0.14920793, -0.14920793, -0.0746039 ], shape=(3145728,))
- zoom :
- 9
- n_side :
- 512
- n_pix :
- 3145728
- nest :
- True
- source_grid_spec :
- HEALPix
Global mean and plot
Computing the global air temperature mean (at the first timestep) and also having a quick plot of it would be a good idea to have as references to compare the upcoming analyses & visualizations to them:
%%time
print(
"Global air temperature average on ", uxda.time[0].values, ": ", uxda.isel(time=0).mean().values, " K"
)
Global air temperature average on 2020-01-02T00:00:00.000000000 : 286.3096 K
CPU times: user 529 ms, sys: 296 ms, total: 825 ms
Wall time: 2.73 s
%%time
projection = ccrs.Robinson()
uxda.isel(time=0).plot(
projection=projection,
cmap="inferno",
features=["borders", "coastline"],
title="Global temperature",
width=700,
)
CPU times: user 16.5 s, sys: 1.13 s, total: 17.6 s
Wall time: 16.4 s
WARNING:param.GeoOverlayPlot00477: Due to internal constraints, when aspect and width/height is set, the bokeh backend uses those values as frame_width/frame_height instead. This ensures the aspect is respected, but means that the plot might be slightly larger than anticipated. Set the frame_width/frame_height explicitly to suppress this warning.
Rasterized point plots
When working with a higher-resolution dataset at a global scale, it’s not always practical to render each cell as a polygon. Instead, we can rasterize the center of each pixel.
projection = ccrs.Robinson()
# Controls the size of each pixel (smaller value leads to larger pixels)
pixel_ratio = 0.5
uxda.isel(time=0).plot.points(
projection=projection,
rasterize=True,
dynamic=False,
width=1000,
height=500,
pixel_ratio=pixel_ratio,
cmap="inferno",
title=f"Global Air Temperature, pixel_ratio={pixel_ratio}",
)
If we decrease the size of each pixel (by setting the pixel ratio to a higher value), we can start to see missing values, which is due to a lower density of points near the poles, leading to some pixels not containing any of our original points.
Because of this, it’s usefully to try a few pixel_ratio values and see which one works best for your given resolution.
projection = ccrs.Robinson()
# Controls the size of each pixel (smaller value leads to larger pixels)
pixel_ratio = 2.0
uxda.isel(time=0).plot.points(
projection=projection,
rasterize=True,
dynamic=False,
width=1000,
height=500,
pixel_ratio=pixel_ratio,
cmap="inferno",
title=f"Global Air Temperature with bad pixel size selection, pixel_ratio={pixel_ratio}",
)